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Abstract 

We improve the current predictions of the standard model for neutron beta decay 
observables and compare them with the currently available ones. Next we study 
their implications in the possible detection of new physics. We discuss where the 
limitations are and where further efforts should be directed. 



Introduction. As precision measurements in neutron beta decay (n/3d) have 
become available, interest in detecting physics beyond the standard model 
(SM) in this decay has increased steadily [1]. However, the predictions of the 
SM for n(3d observables seem to be severely afflicted by our current inability to 
calculate the leading form factors ratio A= g\/ f\ and the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element V u a- It imposes on them two general restric- 
tions, namely, the V — A structure of the weak vertex and the unitary of the 
CKM matrix , respectively. Otherwise these two quantities remain as free pa- 
rameters to be determined from experiment. With current experimental error 
bars a region in the (A, R)- plane (R is the decay rate) can be determined 
where the SM remains valid at, say, 90% CL. We shall refer to this region 
as the standard model region (SMR). Our main purpose in this paper is to 
reduce this region as much as possible, thereby, improving the predictions of 
the SM. Next we shall study the usefulness of this improved region in either 
detecting new physics or imposing more strict bounds on its existence. Finally, 
we shall discuss how future efforts in this respect should be directed. Another 
important aspect to discuss is the limitations for the detection of new physics 
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in nf3d. As precision measurements improve, the limitations will be come sig- 
nificant enough as to render the detection of new physics in this decay very 
obscure or even hopeless. 

Procedure. As we just mentioned agreement of the SM with n/3d can occur 
only in the region of the (A, V u d) plane or equivalently of the (A, R) plane 
where both the V — A and the unitarity restrictions are satisfied, namely the 
SMR. Disagreement, and along with it the possibility to detect new physics, 
corresponds to the region in such planes where the restrictions are no longer 
satisfied. We shall refer to this other region as the new physics region (NPR). 
we shall proceed determine these two regions. 

One can show that up to a precision of 10~ 4 the observables in nf3d depend 
only on two parameters, namely, the (CKM) matrix element V u d and the ratio 
A [2] . All other contributions up to this precision are well determined (through 
the use of the conserved vector-current hypothesis, model-independent and re- 
liable estimates of model-dependent radiative corrections, etc.). We shall use 
the angular asymmetry coeficients a e , a ei/ , and a v (imposing the V — A re- 
striction through A), the CKM V us and V u b matrix elements (imposing the 
unitarity restriction V uh = \Jl — V^ d — V^ s ), and R. Then we form a x 2 func- 
tion with six terms. The first four terms compare the theoretical expressions 
of the a e , a eu , a u and R as functions of A and V u d with their experimental 
counterparts, the fifth term compares V u b of the unitarity square root with its 
experimental counterpart, and the last term constrains V us to its experimental 
value. 

The precision measurements in n/3d occur in R and a e . The state-of-the-art 
will allow that their error bars be reduced to about one tenth of their current 
values [3]. In order to determine the SMR we shall concentrate on these two 
quantities. We shall keep the other four fixed at their present experimental 
values [4], namely, a ev = -0.0766(36), a„ = 0.9830(40), V uh = 0.0036(10) and 
Vus — 0.2196(23). We have folded the theoretical error into the current value of 
R to give R = 1.12879(110) [2]. It is customary to quote the experimental value 
of A instead of a e and we shall follow this practice, i.e., in x 2 we shall replace a e 
by A. There are at present [5,6,7,8] four precise measurements of A (all them 
through a e ), namely, X L = 1.2660(40), X B = 1.2620(50), Ay = 1.2594(38), 
and Xr = 1.2735(21). This last value is slightly revised in [1]. The first three 
are statistically compatible and their average is Xlyb = 1.2624(24). One may 
quote the average of the four, Xlybr = 1.2687(16). However, this average is 
not statistically consistent and one should really keep Xr separate from the 
other three. At any rate the error bars of Xlyb and Xr are comparable, so one 
may take a typical <j\ = 0.0024. Because of this situation we shall cover 5a x 
to the right and 3a to the left of Xlyb- 

For definiteness, we shall work in the (A, R) plane. Our procedure to deter- 
mine the SMR is to calculate x 2 & t points (A, R) covering the rectangle in this 
plane with sides (1.256, 1.274) for A and (1.12549, 1.13209) for R. This latter 
is the ±3<tr range. The values of A and R in the points (A, R) are used as 
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the "experimental" central values in the x 2 function. For the "experimental" 
error bars we shall consider three possibilities. First we take o\ and or at their 
current values, next we take for them one tenth of those current values, and 
last we take for them one thousandth of such values. 

Results. Our numerical results are displayed in Table 1. There we show the 
90% and 95% ranges for A at sample values of R (which go in steps of one 
current <jr up to ±3o\r around its present experimental central value). For 
each value of R we give three entries; the upper, middle, and lower ones cor- 
respond to current error bars, future one tenth error bars, and, so to speak, 
"ideal" one thousandth error bars, respectively. All these numbers are easily 
visualized in Figs. 1-3. The bands within the dotted lines correspond to the 
90% SMR in the three cases considered. In each of these figures we have also 
depicted with I, II, and III the 90%CL ellipses around Xlyb, Xr and Xlybr, 
respectively. In figures 2 and 3 we also show the 90%CL ellipses corresponding 
to one tenth of cr A and a R around these three central values and also around a 
central value previously reported in [9], Xw = 1.2657(30). These figures make 
it very easy to reach several conclusions. The Xlyb and Xr with their current 
error bars both fall completely at a 90%CL in the NPR, i.e., each one gives 
a strong signal for the existence of new physics. This cannot be seen in Fig. 
1, but it becomes very clear in Figs. 2 and 3. Of course, since they fall on 
opposite sides of SMR, the situation is yet undecided. Another conclusion is 
that the SMR is greatly improved by reducing <7r and <j\ by one tenth; but, 
even though there is some improvement, to further reduce <jr and <7\ does not 
provide a substantial gain. The reason is that in Fig. 2 The error bars of V us 
already dominate. This is clearly seen in Fig. 3, where the width of the SMR 
is due almost entirely to the V us error bars. 

The the detection of new physics in nf3d. The ability of n(3d to resolve 
new physics depends on where in the (A, R) plane the future (\ e xp., Rexp.) 
will lie as well as on their future accompanying error bars. One may appre- 
ciate this by looking at the values of \ 2 at points covering the rectangle we 
are studying. For this purpose we have produced Table 2. In each entry we 
give two numbers. The upper one uses the current or and <J\, the lower one 
uses one tenth of or and <j\. \ 2 does not further increases substantially if one 
uses one thousandth of a R and <J\, but we shall not produce the corresponding 
numbers in Table 2. x 2 ma y reach impressively high values. 
There is another way to appreciate this resolving power. This can be seen in 
Table 3. We selected the four values of A marked as I-IV in Figs. 2 and 3. The 
entries in this table contain four values of \ 2 ■, which up-downwards are in the 
order I-IV. The error bars in A and R are reduced from their present values 
to one tenth of them in steps 1/N, N = 1, 2, 10. Each one of the A values 
corresponds to new physics, as is indicated by the small ellipses in Figs. 2 and 
3. In Table 3 we can see how the resolution increases as the error bars in A 
and R decreased. 

To clearly resolve I and II it is quite sufficient to cut or and a\ simultaneously 
to a little over one-half their current values. To resolve point IV requires that 
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or and <7\ be simultaneously cut to around one sixth. To resolve the overall 
average of A at III requires more effort , because of its closeness to the SMR. 
However, Fig. 3 shows it falls in the NPR. Also, one can see in this table that 
reducing a R alone is of little use in detecting new physics. While reducing a x 
alone helps in cases I and II, but is useless in cases III and IV. 
Summary and discussion. In this paper we have substantially improved 
the predictions of the SM for nj3d observables, as can be appreciated by com- 
paring Fig. 1 with Figs. 2 and 3. Our analysis has several implications. First, 
it allows to clearly and also easily appreciate recent precision measurements, 
XLYB and \ R . In Fig. 1 they show some indication in favor of new physics, 
while in Figs. 2 and 3 each one is clearly resolved past the 90% C.L. in favor 
of the latter. Of course, we still have to wait for a consistent average to be 
determined in the future . Second, our current inability to compute A and R 
with the SM is not a real obstacle to detect new physics in nf3d. Ideally, the 
SM predictions should be a point in the (A, R)- plane. Instead, with A and 
V u d remaining as free parameters, the predictions become still ideally a line 
in this plane. Thus, it is still possible to resolve this line from new physics 
contributions. Third, the real limitations for this enterprise come from the 
width of the band around this ideal line. As we have seen, there are only two 
important ingredients in this width. One is the 10" 4 theoretical uncertainty in 
R and the other one is the error bars on V us . The first one turns out to be not 
important. It forces us to stop at the band of Fig. 2. However, looking at Fig. 
3, which assumes that this theoretical error has been reduced to zero, we learn 
that error bars of V us are the real limitation for nf3d to detect new physics. 
In the light of this discussion, we can understand better where future efforts 
should be directed. On the theoretical side, they should be directed towards 
reducing the theoretical uncertanties in R and in the determination of V us . On 
the experimental side, our analysis makes it clear (see Table 3) that precision 
measurements must simultaneously reduce the error bars on R and A. Improv- 
ing measurements in only one of them is useless. The precision measurements 
of the other two angular coefficients a eu and a u lead to the same bands of 
Figs. 2 and 3. In this sense they are no longer needed. However, if new physics 
is indeed detected in nf5d, they will be very useful in delucidating what kind 
of new physics it would be. Finally, it is clear that the determination of V us 
must be greatly improved, because it represents the real limitation for nf3d to 
detect new physics. 

Acknowledgments. The authors are grateful to CONACyT (Mexico) for 
partial support. 
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Table 1 

Ranges for A at fixed values of R, in the three cases for the size of cr\ and ctr 
discussed in the text. For completeness, we also include the range for V u a- Xo * s the 
minimum of \ 2 an d ^0 is the corresponding A. 



Rn 



xj 



Ao 



bounds 
C.L. 



o n 



A 



bounds on V U( i 
C.L. 90% C.L. 



(1.26581, 
1.26911, 
1.26952, 

1.26503, 
1.26835, 
1.26877, 

(1.26424, 
(1.26761, 
(1.26802, 

(1.26345, 
(1.26687, 
(1.26727, 

(1.26268, 
(1.26611, 
(1.26651, 

(1.26189, 
1.26536, 
1.26576, 

(1.26110, 
(1.26461, 
(1.26501, 



1.27464) 
1.27072) 
1.27041) 

1.27386) 
1.26997) 
1.26966) 

1.27308) 
1.26922) 
1.26891) 

1.27231) 
1.26846) 
1.26816) 

1.27154) 
1.26770) 
1.26741) 

1.27076) 
1.26695) 
1.26666) 

1.26998) 
1.26620) 
1.26591) 



(1.26471 
(1.26883 
(1.26940 

(1.26393 
1.26807. 
1.26865, 



1.27574) 
1.27103) 
,1.27060) 

1.27496) 
1.27029) 
1.26985) 

1.27418) 
1.26953) 
1.26910) 

1.27341) 
1.26876) 
1.26835) 



0.97270) 
0.97519) 
0.97597) 

0.97272) 
0.97519) 
0.97597) 

0.97274) 
0.97519) 
0.97597) 

0.97276) 
0.97520) 
0.97597) 

0.97277) 
0.97521) 
0.97597) 

0.97279) 
0.97521) 
0.97597) 

0.97281) 
0.97521) 
0.97597) 



1.13209 



1.13099 



1.12989 



1.12879 



1.12769 



1.12659 



1.12549 



0.71 
0.75 
0.75 

0.65 
0.71 
0.69 

0.61 
0.67 
0.64 

0.56 
0.64 
0.59 

0.54 
0.60 
0.56 

0.52 
0.58 
0.52 

0.50 
0.56 
0.50 



1.26988 
1.26988 
1.26990 

1.26913 
1.26913 
1.26915 

1.26839 
1.26839 
1.26840 

1.26763 
1.26763 
1.26765 

1.26689 
1.26689 
1.26689 

1.26614 
1.26614 
1.26614 

1.26539 
1.26539 
1.26539 



(1.26314. 
(1.26734. 
(1.26790. 

(1.26235. 
(1.26659. 
(1.26715. 



(1.26158. 
(1.26582. 
(1.26639, 

(1.26079, 
1.26508, 
1.26564, 

(1.26001, 
(1.26433, 
(1.26489, 



1.27264) 
1.26801) 
1.26760) 

1.27186) 
1.26725) 
1.26685) 

1.27108) 
1.26652) 
1.26610) 



(0.97832. 
(0.97621. 
(0.97539, 

(0.97834. 
(0.97622, 
(0.97539, 

(0.97837. 
(0.97622, 
(0.97539, 

(0.97840. 
(0.97622, 
(0.97539. 

(0.97842 
(0.97622. 
(0.97539. 

(0.97845. 
(0.97623. 
(0.97539, 

(0.97848. 
(0.97623. 
(0.97539. 
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Table 2 

Values of x 2 at sample points in the (A, R) plane for current o\ and <tr (upper 
entries) and for one tenth of their current values (lower entries) . 



JX \A 


1 9^9 


1 9^7fi 


1 9^nn 

l.ZOUU 


1 9^9/L 


1 9fi/IR 


1 9^79 


i.zoyo 


1 9790 
1 .Z / zu 


1 97/M 


1.13209 


33.7 


23.7 


15.7 


9.2 


4.7 


1.9 


0.8 


1.3 


3.4 




3352 


2353 


1440 


805 


354 


87 


1.1 


32 


230 


1.13099 


30.4 


21.0 


13.4 


7.6 


3.6 


1.3 


0.7 


1.7 


4.4 




2924 


1983 


1224 


647 


253 


41 


0.8 


74 


328 


1 12989 


27 3 


18 4 


11 4 


6 2 


2 7 


9 


0.8 


2 3 


5 6 




2613 


1728 


1026 


506 


168 


12 


5 


133 


443 


1.12879 


24.3 


16.0 


9.5 


4.9 


1.9 


0.7 


1.0 


3.1 


6.8 




2320 


1491 


845 


381 


100 


1.2 


24 


208 


675 


1.12769 


21.5 


13.8 


7.8 


3.8 


1.3 


0.6 


1.5 


4.0 


8.2 




2044 


1272 


682 


275 


51 


0.6 


62 


302 


725 


1.12659 


18.9 


11.7 


6.4 


2.8 


0.9 


0.6 


2.0 


5.1 


9.9 




1784 


1069 


536 


185 


17 


2.4 


117 


414 


893 


1.12549 


16.4 


9.9 


5.1 


2.0 


0.6 


0.9 


2.8 


6.3 


11.7 




1543 


884 


407 


114 


2.1 


18 


190 


543 


1079 
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Table 3 

Increase of x 2 as a\ and or are reduced as 1/N, N = 1,2, 10, at the points I- 
of figures 2 and 3. The up-downwards order is also I— IV. 



IV 



R\X 


H 

cr\/l 


oa/2 


/o 

oa/3 


/ A 

crx/4: 


oa/5 


oa/6 


in 

o-a/7 


/o 

oa/8 


oa/9 


/-i r\ 

ax/W 


a R /l 


4.8 


13.1 


21.1 


27.3 


31.6 


34.6 


36.7 


38.5 


39.3 


40.1 




7.4 


20.2 


31.5 


39.3 


44.2 


47.6 


50.0 


52.0 


52.9 


53.8 




0.8 


1.0 


1.1 


1.1 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 




1.2 


2.3 


3.1 


3.8 


4.3 


4.6 


4.8 


5.0 


5.1 


5.2 


a R /2 


5.1 


16.3 


31.9 


49.3 


66.2 


81.5 


94.9 


106 


116 


123 




8.0 


26.6 


51.2 


77.0 


101 


121 


139 


152 


164 


173 




0.8 


1.1 


1.2 


1.3 


1.5 


1.5 


1.6 


1.6 


1.7 


1.7 




1.2 


2.6 


4.2 


5.9 


7.7 


9.4 


10.8 


12.0 


13.1 


14.0 




5.2 


17.1 


35.3 


58.1 


83.4 


109 


135 


160 


182 


203 




8.1 


28.1 


58.2 


94.5 


133 


172 


209 


242 


272 


298 




0.8 


1.1 


1.3 


1.5 


1.6 


1.7 


2.0 


2.1 


2.2 


2.3 




1.3 


2.6 


4.5 


6.8 


9.4 


12.1 


14.8 


17.3 


19.8 


22.0 


an/ A 


5.2 


17.5 


36.7 


62.5 


91.8 


125 


159 


194 


228 


262 




8.1 


28.8 


61.2 


103 


150 


202 


254 


305 


354 


400 




0.8 


1.1 


1.3 


1.5 


1.7 


1.9 


2.1 


2.3 


2.5 


2.6 




1.3 


2.7 


4.7 


7.1 


10.1 


13.5 


17.0 


20.6 


24.9 


27.8 




5.2 


17.6 


37.7 


63.9 


96.2 


133 


173 


215 


259 


303 




8.1 


29.1 


62.7 


107 


160 


219 


282 


347 


412 


476 




0.8 


1.1 


1.3 


1.5 


1.8 


2.0 


2.3 


2.6 


2.8 


3.1 




1.3 


2.7 


4.8 


7.5 


10.8 


14.6 


18.3 


22.7 


27.8 


31.7 




5.2 


17.6 


38.0 


65.1 


98.8 


138 


182 


230 


279 


330 




8.1 


29.3 


63.5 


109 


166 


230 


300 


375 


452 


530 




0.8 


1.1 


1.3 


1.5 


1.8 


2.1 


2.4 


2.8 


3.1 


3.4 




1.3 


2.7 


4.7 


7.6 


10.8 


14.7 


19.2 


24.6 


29.8 


34.3 


a R /7 


5.3 


17.7 


38.3 


65.8 


101 


141 


189 


238 


293 


350 




8.1 


29.4 


64.0 


111 


169 


240 


312 


394 


480 


569 




0.8 


1.1 


1.3 


1.6 


1.9 


2.2 


2.5 


2.9 


3.3 


3.7 




1.3 


2.7 


7.7 


11.2 


11.2 


15.0 


19.7 


25.4 


30.3 


36.1 




5.3 


17.7 


38.4 


66.2 


102 


144 


192 


245 


302 


364 




8.1 


29.4 


64.4 


112 


172 


242 


321 


408 


501 


598 




0.8 


1.1 


1.3 


1.6 


1.9 


2.2 


2.6 


3.0 


3.4 


3.9 




1.3 


2.7 


4.7 


7.7 


11.3 


15.2 


20.0 


25.4 


32.0 


37.4 




5.3 


17.7 


38.5 


66.6 


102 


145 


194 


249 


309 


374 




8.1 


29.5 


64.6 


113 


173 


245 


327 


418 


516 


620 




0.8 


1.1 


1.3 


1.6 


1.9 


2.3 


2.7 


3.1 


3.6 


4.0 




1.3 


2.7 


4.7 


7.7 


11.1 


15.4 


20.3 


25.8 


31.9 


39.3 


VR./W 


5.2 


17.8 


38.3 


66.8 


103 


146 


196 


253 


315 


381 




8.1 


29.5 


64.8 


113 


175 


248 


332 


425 


527 


636 




0.8 


1.1 


1.3 


1.6 


1.9 


2.3 


2.7 


3.2 


3.6 


4.2 




1.3 


2.7 


4.8 




11.2 


15.5 


20.5 


26.1 


33.2 


39.1 
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R 



1.132 F~r 
1.131 - 
1.13 
1.129 
1.128 
1.127 - 
1.126 - 



_L 




1.256 1.258 1.26 1.262 1.264 1.266 1.268 1.27 1.272 1.274 

A 

Fig. 2. The SMR, at one tenth of the current values of a\ and ctr, is the band within 
the dotted lines. See text for other explanations. 




Fig. 3. The SMR, at one thousandth of the current values of a\ and or, is the band 
with the dotted lines. See text for other explanations. 
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